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Abstract 

We evaluate the two-photon box contribution to heavy fermion production in electron 
positron annihilation, that provides 0(a 2 ) electromagnetic corrections to the Born 
cross section. The study of its non-relativistic expansion, relevant at energies close 
to the threshold of production, is also performed. We also verify that the threshold 
expansion of the one-loop integrals correctly reproduces our results, thus extending the 
applicability of this technique to heavy fermion production diagrams. 
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1 Introduction 



Heavy fermion production processes out of electron positron annihilation, e + e~ — > //, 
have become a subject of thorough study in the last years. Their interest embodies multiple 
features and a wide energy range, from close to threshold production to high-energy colliders. 
LEP and LEP2 have provided the appropriate tool pushing behind this burst. In addition 
this is among the scattering processes with higher expected number of events at a future 
Linear Collider running in the 0.5 TeV — 1 TeV energy region like TESLA and NLC/JLC-X, 
or CLIC at higher energies. Their interest arises mainly from the possibility of exploring 
New Physics and, therefore, an accurate description within the Standard Model is necessary 
for the analyses of data. Projects like ZFITTER |L| and the ongoing CalcPHEP [^, |3| aim 
to provide the relevant theoretical framework for that purpose. 

QED corrections seem to be of little interest when probing the quantum effects within 
the Standard Model, but it is obvious that their contribution, however small, should be 
considered in order to disentangle New Physics effects. Besides, if a deeper understanding 
on the physical parameters of heavy fermions is intended, electromagnetic r + r~ and heavy 
quark QQ production out of e + e~ annihilation at threshold energies supplies the required 
information. 

From a theoretical point of view e + e~ — > ff cross sections close to threshold evaluated 
within perturbation theory are mislead due to the presence, in the physical system, of a 
kinematical variable of the same order than the gauge theory coupling : the velocity of the 
heavy fermion pair in the center of mass of the colliding system, (3 = J 1 — 4M 2 /s, with M 
the mass of the / fermion. Hence, when (3 <~ a, care has to be taken in order to resummate 
terms as (a//3) n or (a\iaj3) n that can give potentially large contributions ||. Recently the 
development of non-relativistic effective field theories of QED and QCD |J implements the 
suitable systematic procedure to follow. Facilities as the proposed Tau-Charm Factory, 
a high-luminosity e + e~ collider with a center-of-mass energy near the r + r _ production 
threshold ||, would provide excellent information on the mass of this lepton 0. Moreover 
an accurate determination of the mass of the top quark (difficult to get at the next hadron 
colliders) requires a future lepton collider at the tt threshold 0. Consequently a thorough 
study of the non-relativistic contribution to cr(e + e~ — > ff) both from electromagnetic and 
strong interactions is mandatory. 

In Ref. a detailed study of the threshold behaviour of cr(e + e~ — > t + t~) was performed, 
and it was pointed out that, within the 0(a 2 ) electromagnetic corrections to the Born cross 
section, the squared amplitude of the box diagram involving two-photon r + r~ production 
(see Fig. |I]) had not been considered yet. The electroweak one-loop contributions to the 
e + e~ — > ff process were evaluated in Ref. ||. Here this box contribution was already taken 
into account, though an explicit expression was only given for the M = case. In this paper 
we provide the amplitude of this diagram for a final massive fermion. [] 

Once the explicit result is worked out we perform its non-relativistic expansion in terms 

1 While writing this article Ref. || appeared. In this preprint a full expression for the QED box diagram 
amplitude is also given. 
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Figure 1: Direct (a) and crossed (b) box diagrams for e + e — > //. 



of the (3 velocity and we find that the contribution of this diagram to the cross section is of 
0(a A f3 3 ), that is 0(a 2 f3 2 ) over the Born cross section. The additional suppression driven by 
the velocity squared indicates that the contribution of the two-photon box diagram to the 
production of heavy fermions at threshold is negligible compared to the precision foreseen in 
the next future. 

In Section 2 we detail the calculation of the box diagram contributing to e + e~ — > // in 
the limit when m e <C M, and we provide the full analytical result. Section 3 is dedicated 
to the study of the threshold behaviour of the box amplitude as obtained directly from 
our previous result. We confirm the features of this threshold amplitude by performing an 
alternative analysis of the integrals through the asymptotic expansion method in Section 4. 
Our conclusions are collected in Section 5. Finally, two Appendices contain the basic scalar 
integrals appearing in the article and a comment on the infrared divergent part of the box 
amplitude. 



2 Two— photon box diagram 

The contribution to the 5-matrix of the process e~(p)e + (p') —>■ f(k)f(k') of the two- 
photon box amplitudes is depicted in Fig. [I] and it is defined by 

(ff\tT\e + e-) box = (2ir)U(p + p'-k-k f )tM box . (1) 

As we are interested in heavy fermion production we will perform the evaluation for 
k 2 = k' 2 = M 2 and p 2 = p' 2 = m 2 <C M 2 (we neglect the electron mass where possible). The 
two-photon box amplitude is gauge invariant and, consequently, we perform the calculation 
by taking the Feynman choice. The direct box amplitude, Fig. 0(a), is written down following 
QED Feynman rules as : 

M =e 4 2 f dH { ^ p,) tlj^PR ^ g ~ t + I + M ) It Mk')} (?) 

^ f J i(2vr) 4 (P - m 2 ) [(£ + k - p) 2 - M 2 ] [(£ - p) 2 - A 2 ] [(£ + p') 2 - A 2 ] ' 1 ' 

where we have introduced a photon mass A in order to regularize the infrared divergences 
known to be present in this amplitude. The crossed box diagram in Fig. ffl(b) can be obtained 
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from (fj) by adding an overall minus sign, reversing the order of the 7^, 7„ matrices in the 
heavy fermion bilinear, and performing the substitutions k — > k! everywhere (except for 
the spinors) and M — ► —M. Hence, in Eq. ([[]), M.box = M-a + -Mb- The evaluation of 
the integrals is slightly cumbersome but straightforward and the details are given in the 
Appendix A. 

With the definition of the Mandelstam variables s = (p+p') 2 and t = (p—k) 2 , the spinorial 
structure of M.box is decomposed into 4 sets of amplitudes Lf K multiplied by corresponding 
coefficients tof : 

M box {^s,t) = e i Q)Y J E L T< , (3) 
i=i p=±i 

with the basic amplitudes 

LT = [v e {p') lfl P K u e {p))[u f (k)r{l + Kp l5 )v f (k')}, 
Lf = [v e (p')fkP K u e (p)] [u f (k)p , (l + Kp l5 )v f (k')} , 
Lf = [*.(?/) *P B ti.(p)] [u f (k)(l + K P ^)v f (k')] : 

LT = [v e (p') 1 ,P K u e (p)}[u f (k) 1 ^(l + K Pl5 )v f (k')} . (4) 
The latter have been written in terms of the initial state e + e~ chiral projectors 

P, = ^(1 + ^75) , k = ±1 , (5) 

which, as we are considering massless initial fermions, satisfy P K u e (p) = u e (p), being k the 
initial electron helicity; in the massless limit, positron helicity is forced to be — k to have 
a non vanishing amplitude. The dependence of Adbox on the spin state of the final state 
fermions has not been explicitly stated. 

The w? coefficients can be written in terms of four auxiliary functions 2 = 0,1,2,3: 

wf = ]-T Q {s,t) , 

w\ = ~ ^(s, u) , 

w t = Fi{s,t) + F-z{s,t) , 

w 2 — Fi{s,u) + F%{s,u) , 

w£ = M(f x {s,u)- Fi(s,t) + F 3 (s, u)-F 3 (s,tj) , (6) 

W3 = m(f z {s,u) - Fz{s,t)) , 

w+ = --M^ 2 (s,t)-^ 2 (s,u)) , 
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that read 



M 2 — t ( M 2 — s — / 

F (s,t) = Ffo,t) + ^ jX9 , J t(sD - 2C t )-M 2 ' 



(M 2_ t) 2 + sU -- U ~w ~ M2 _ t (*A>-2C t ) (7) 
_ (2M 2 - s - 2t) (C s + C M ) } , 



4M 2 
M 2 — t 



{B M - B t ) 



(8) 



+ 



M 2 -t 



(M 2 - t) 2 + st 



(2M 2 -s-2t) ((M 2 - t) (sD Q - 2C t ) + s C s 
+ ( 2(M 2 - tf + s(2t + s - 4M 2 ) ) C A/ } , 



F 2 {s,t) = Ffo,t) - (M2 _ 1 t)2 + gt {(^ 2 + ^) (sD -2C t )- {s-AM 2 )C 



M 



(9) 



F 3 (s,t) 



1 



(M 2 - t) 2 + I M 2 - t 

(M 2 - 1) 2 



- (2M 2 - s - 2t) C S J 



2t fa » ^ + 2 ( M2 + t ) r p r ^ 

2 — 7 \£>t — d M ) H — 7^ — — \a s — ±s Mj 



+ 



(M 2 - t) 2 + st 
1 



am* - s 



(M 2 -t) ( Sj D -2C t ) + sC s 



(10) 



(4M 2 -s)((M 2 -t) 2 + si) 



-4M 8 + 6(s + 2t)M 6 

- (s + 2t) (s + 6t) M 4 

■f 2t (s 2 + ts + 2t 2 ) M 2 + s 2 * 2 



C 



Full expressions for the scalar functions B s , B t , Bm, C s , C t , C t , Cm, A) and D can be found 
in Eqs. ( |A.5[ - |A.1^D of the Appendix A. It can be seen from ^(s, t) and ^(s, t) that the two- 
point functions B t , B s and 5^ only appear in non-divergent combinations while the rest 
of scalar integrals in ^(s, t) are UV finite. Clearly M.box{^\ s,t) is ultraviolet finite. Scalar 
integrals D (-Do) and Ct (Ct) are infrared divergent for vanishing photon mass A, however 
the combinations sDq — 2Ct and sDq — 2Ct are divergenceless. Hence all the divergences in 
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the !Fi(s,t) functions are collected in ^(s, t) and J 7 ^ (s,t), given by : 

F£( 8 ,t) = 2 d , (ii) 

and we get for the infrared divergent part of the box amplitude : 

A more complete discussion on the infrared structure of the QED box diagram and the 
determination of is relegated to Appendix B. 

Incidentally our result can be used to evaluate the similar two-gluon box contribution 
to the heavy quark production out of light quarks, q(A)q(B) — > Q{C)Q(D) (between paren- 
theses we label the colour quantum numbers). In order to get this amplitude we need to 
substitute the e A Qj factor in Eq. (|) by gj(t b t a ) BA {t a , t b } CD (a sum over repeated indices is 
implied). 

We have checked that our amplitude in Eq. (^), when summed over polarizations, coincides 
with a recent result found in Ref. ||, though these authors use a different basis of spinor 
operators. Moreover, from our calculation for Aibox, we can recover the case where the final 
fermions are massless. The limit M — > can be directly applied to the wf coefficients, 
Eqs. (@), and to the scalar integrals quoted in the Appendix A. Within this limit our result 
agrees with the earlier calculation in Ref. ||. 

3 Heavy fermion production at threshold 

Close to // threshold, it is more convenient to expand the production amplitude in terms 
of the fermions velocity in the center-of-mass of the colliding system (3 = Jl — 4M 2 /s. 
Hence production amplitudes are written in a combined expansion in powers of a and f3, 
and the importance of each contribution is estimated taking a ~ (3. This feature spoils the 
perturbative expansion in QED due to the appearance of 0(a n /(3 n ) and 0(a m ln n f3) terms 
that diverge as (3 — > 0. As a consequence, a resummation of such terms is necessary to avoid a 
breakdown of the perturbative series, and well-known results from the familiar non-relativistic 
quantum mechanics are obtained. Nevertheless it is somewhat misleading to associate the 
appearance of these Coulomb terms to the non-relativistic motion of the fermion pair, as 
the scattering amplitude calculated from quantum mechanics does not show any kinematic 
singularity close to threshold: their ultimate origin is the inadequacy of the diagrammatic 
QED expansion in powers of a to account for the correct non-relativistic dynamics. Keeping 
this in mind, one should not discard, a priori, divergent terms in the velocity appearing in 
any QED diagram involving fermions with small velocities. 

2 In this colour factor t* = A l /2, where A* are the SU(3) Gell-Mann matrices and Tr(t' l t j ) = ^ jl. 
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In Ref. [[F) it was pointed out that the contribution at threshold of the two-photon box 
diagram should be analysed in a NNLO calculation of cr(e + e~ — > t + t~). In this Section 
we proceed to perform the expansion on M.b ox as given in Eq. (|3|). The leading terms in 
the velocity expansion of the coefficients wf can be obtained by taking into account the 
dependence of the Mandelstam invariants s, t, u on the velocity (3 and the angle 9 between 
the momenta of the heavy fermion and the electron in the colliding center-of-mass system. 
The relation is given by : 



4M 2 
1 -(3 2 



t = M 2 



2M 2 
l-(3 2 



[I — (3 cost 



u = M 2 



2M 2 
l-(3 2 



[l + f3cos9) . (13) 



Carrying these expressions to the wf coefficients displayed in Eq. @ and neglecting 0{(3 2 



terms we obtain : 
1 

w 



1 "~ 384M 2 tt 2 



ix 2 + 3 In 



2 4M 2 2 m 2 
A 2 "X 2 " 



8- 14m -8 In 2 + 12 In 



AM 2 



(3 cos 9 



+ 0{(3 2 ) , 



w n 



-«,+(/?--/?) 
1 



384M 4 vr 2 



7T -8 + 8zvr + 81n2-31n^ 



4M 2 , 9 m 2 

— +31n V 



+ ^r 2 -34 + 4ivr + 161n2 + 121n-^- 

+ 0(f3 2 ) , 



3\n 2 



AM 2 
"A 2 " 



2 

+ 31n 2 — )(3cos9 



Wn 



wl((3 - -/3) 



(14) 



240M% 2 
-1 



37 + 2i7r-641n2 )(3cos9 + 0{/T) 



w 



'• "' 480M% 2 
-1 



ll + OT-321n2 I3cos9 + O(0 



w 



• " 384M% 2 



7r 2 + 627r-48 ln2 + 121n 



4M 2 



3 In 2 



4M 2 



m 



X 2 X 2 



+ 3\n z — )f3cos9 + 0((3 2 ) , 



w , 



-w A 



The amplitudes Lf K , containing fermion wave functions, must also be expanded in terms of 
(3 to fulfill the expansion of M.box at small velocities. We shall not give the full result of such 
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expansion, but just quote their leading behaviour, which can be easily obtained by choosing 
an explicit representation of the gamma matrices and spinors. We thus get: 



Lf = 0(1) 



Lf = 0((3) , 



Lf = 0(1) . 



(15) 



The terms quoted in Eqs. ( |14| ) together with the expansion in Eq. (|T5| ) allow us to obtain the 
leading near threshold contribution to the cross section of the box amplitude Aibox- Recall 
that, by virtue of Furry's theorem, the interference of the QED box amplitude with other 
one-loop amplitudes for the process e + e~ — > // vanishes and, consequently, 1 f, 0;E | 2 adds 
incoherently to the rest of 0(a 4 ) corrections to <x(e + e~ — > //), as studied in Ref. fTj. The 
final result for the squared and averaged box amplitude is : 



1 



box 



(Qf 



a) 



pol. 




In2) : 



(16) 



+ 



J M 



24-12< m + ?ln2 + ^ 



M 



1 6 

+ (- 8C+(— ln2 + 



-7T 



320 



2881n2 + -7r^ + 32]L M 



+ 564+(- 2881n2 + V + 32)^ + ^ln 2 2-^ln2 
^ 3 ' 9 9 



7T^_8 

18 9 



14 



7r 2 ln2 7T- 



16 



cos 2 ^ \(3 2 + C(/3 3 ) , 



with 



L M = In 



4M 2 



and 



m 



(17) 



Hence we conclude that the result in Eq. ([IBD , proportional to a 4 j3 2 , represents a N 4 LO 
correction with respect the LO result (the tree level e + e~ — > // amplitude squared, which 
is already of 0(a 2 )). In Ref. [0, box amplitudes were not included with the rest of the 
one-loop diagrams to complete the NNLO calculation of a(e + e~ — > t + t~) at threshold, their 
behaviour with (3 being unknown. Our evaluation of |A^b OI | 2 has proven that this is, indeed, 
j3 2 suppressed with respect the NNLO contributions considered in M. 



4 Threshold amplitude by asymptotic expansion of in- 
tegrals 

The counting of powers of the velocity appearing in a defined amplitude is not straight 
because (5 is not a parameter in the Lagrangian, but rather a dynamic scale which is driven 
by the propagators inside loop integrals. In recent years, this issue made awkward to define 
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a non-relativistic effective theory suitable for describing quarks and leptons at low velocities. 
Important progress was made after the development of the threshold expansion by Beneke 
and Smirnov JTIJ. This technique allows for an asymptotic expansion of Feynman integrals 
near threshold, providing a set of much simpler integrals which are manifestly homogeneous 
in the expansion parameter and so have a definite power counting in the velocity. The 
procedure should confirm that the two-photon box amplitude is not enhanced at low (3, as 
we have found by explicit evaluation. This we discuss in the following. 

The expansion method, described in Ref. [pTOfj , begins by identifying the relevant momen- 
tum regions in the loop integrals, which follow from the singularity structure of the Feynman 
propagators dictated by the relevant scales that appear in the problem. For on-shell scatter- 
ing amplitudes of heavy fermions, three scales are identified: the heavy fermions mass, M, 
their relative 3-momentum, \k\ ~ M/3 and their energy k Q ~ M/3 2 . Accordingly, the loop 
four momentum near the singularities can be in any of the following regimes: 

hard : £ ~ \£\ ~ M, 

soft : £ Q ~ \£\ ~ M/3, 

potential : £ ~ M/3 2 , \£\ ~ M/3, 

ultrasoft : £ ~ \£\ ~ M/3 2 . (18) 

The original integral is then decomposed into a set of integrals, one for every region, and 
a Taylor expansion in the parameters, which are small in each regime, is performed. Every 
integral, containing just one scale, will thus contribute only to a single power in the velocity 
expansion. The procedure requires the use of dimensional regularization in handling the 
integrals, even if they are finite, in order to assure that the result from each regime just 
picks up the corresponding pole contribution and vanishes outside. Following this heuristic 
rules, the authors of Ref. [TO reproduce the exact (3 expansion of some one-loop and two- 
loop examples. Although a formal proof of the validity of the asymptotic expansion close 
to threshold has not been given, the perfect agreement in the examples supports their use 
in general one-loop diagrams. We provide a new test by addressing the rules to the QED 
box amplitude with e + e~ in the initial state, extending the use of the threshold expansion 
to diagrams with heavy and massless fermions in the external legs (i.e. production-like 
diagrams). We will keep the electron mass finite along the procedure, although much smaller 
than any other scale, to keep track of the logarithms of m present in the box amplitude. 

Our amplitude M.^ ox is characterized, as shown in the Appendix A, by the four point 
integrals D , D^, D^ u in ( |A.1| ). If present, inverse powers of the velocity in Aibox can only 



originate from these integrals. In addition, we can focus on the behaviour of the scalar 
integral D Q , as the £^£^£ v vectors in and will produce factors of one of the scales 
of the problem (M, M/3 or M/3 2 ) in the numerator of the amplitude without affecting the 
leading singular behaviour in (3. Let us change the routing of momenta in D (|A.l ) in order 
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to make the scaling arguments more transparent: 
D, 



*(2vr) D [(Q/2 + T/2 - £) 2 - m 2 } [(Q/2 + R/2 - I) 2 - M 2 ] [I 2 - A 2 ] [(Q - i) 2 - A 2 ] ' 

(19) 

where the standard +i5 prescriptions are implicitly understood in the propagators, the Q and 
R vectors are defined in relation with Eq. ( |A.1| ) and T = p — p'. The external four vectors 
Q and R scale as M and Mf3 respectively, while T 2 = —s + 4m 2 ~ M 2 . Using momentum 
T is preferred to the electron (positron) momentum p (p') because, the spatial and time 
components of the latter, although scale as M, cancel out in the total momentum squared 
p 2 = m 2 ~ 0. The infrared regularization of the integrals is automatically guaranteed by 
dimensional regularization and, therefore, we will not longer retain a fictitious mass for the 
photon. 

In the potential region £ <C \£\ <C M and, accordingly, we can expand terms in the 
propagators. The leading contribution is 

r d D £ 1 
nP = / _ (of)) 

J i(2n) D (£ ■ T) (-f + £ R - Q £ ) (-£ 2 ) (Q§) ' 1 } 

where we have also dropped the term — £ 2 in the electron propagator to be compared to 
£ ■ T ~ M 2 f3. The overall scaling of the potential integration is easily estimated to be of 
order M 4 /3 5 /M 8 /3 5 ~ 1/M 4 , so no velocity enhancement is this region is expected. In fact, 
the integral above is zero because, closing the C.q integration contour in the lower half-plane, 
the pole at £ = (£-R—£ 2 )/Q +iS lies outside Similarly, subleading terms in the expansion 
of propagators in this region are vanishing, as they share the same pole structure. 

When the loop momentum I is soft or ultrasoft, the assumption £q ~ \£\ <§C M leads to 
the same expansion of the propagators in Dq: 



d D i 1 



n s > us = / : (21) 

J z(2tt) d (£-T- Q £o) (-Q £ ) (I 2 - £ 2 ) (Q 2 ) ' 1 ; 

It scales as 1/M 4 in both the soft and ultrasoft regimes and, indeed, vanishes in dimensional 
regularization because, after picking up the residue in the lower plane, £o = \£\ — i5, the 
remaining D — 1 dimension integral is scaleless: 

1 r d D -H 1 1 



2Q 3 J (2tt) d - 1 \£\ 2 (Q \£\-T-£) 
dQp-i 1 



2QIJ (27T> D - 1 (Q - |T| cosyp) 



3 Notice that the ^-integration in Dq does not vanish in the outer semicircle. Rigorously we should 
keep the £q term in the heavy fermion propagator, so Dq is well defined. Poles would then be located at 

£q = I ^<5o i \J Qo ~ 4(-^ • R — f 2 ) — iS^j . The root £q scales as M and is taken into account in the hard 

region while £q = (£• R- 1 2 ) /Q + iS once we consider that \£\ < M in the potential region, and we recover 
the above result. 
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with (p the angle between the vectors T and t. The same argument holds for subleading 
terms in this region. 

Finally, the integral in the hard region is obtained by dropping out terms involving non- 
relativistic fermion three-momenta from propagators. Hence, the only scale which remains is 
the hard parameter M, and so there is no additional velocity dependence in the denominators. 
More explicitly, the expanded integral in the hard regime, at leading order in /3, is 



h,0(l) 




d D i 



i(2ir) D (£ 2 - £ ■ T - Q ■ £) (£ 2 - Q ■ £) £ 2 (Q 



(23) 



and there is no need to separate time from spatial components in the integration. The above 
integral trivially scales a 1/M 4 , and its explicit calculation in D = 4 — 2e dimensions has 
been performed following Ref. 



D 



h,0(l) 



-2e 



In 



In 



+ ln(47r) - 7^ 



(24) 



where je is the Euler-Mascheroni constant. Terms proportional to the electron mass m have 
been dropped. The pole in Eq. ( pll ) is of infrared origin, and it is the analogous to the In A 2 
term in the full result of D , Eq. (|A.5|) . Indeed, Eq. (24) reproduces the leading term in the 
velocity expansion of D , after the usual replacement In A 2 — > (47r) e /T(l — e)/e. 

The following order in the expansion within the hard region would have a £ ■ R = —£ ■ R 
term in the numerator, and it would behave as (3/M 4 : 



D 



h,0(f3) 




d 



R 



z(2tt) u (£ 2 -£-T-Q-£)(£ 2 -Q-£) 2 P (Q - 



R-T 

(3 cos 9 
8n 2 s 2 



D 



d 



i(2n) D (£ 2 -Q- 



fx 



-2c 



In 



In 



) 2 £ 2 (Q- 
—s — i5\ 



+ ln(47r) - 7 £ 



(25) 



which agrees with the second term in the velocity expansion of Dq. The series expansion in 
(3 of the scalar function Dq is thus reproduced by that of Dq, while the rest of integration 
regions does not contribute at all. 

Therefore we have seen, by asymptotic expanding the integral before its computation, that 
the box amplitude receives no contributions from the regions of potential, soft and ultrasoft 
loop momentum, and it is then preserved from Coulomb type singularities, as it was shown 
by explicit calculation. This fact reveals that, as expected, the box production graph is a 
process dominated by the high scale, as it involves annihilating photons which carry energies 
of the order of the mass of the non-relativistic fermions. 

Let us finally note that, although we have reproduced the (logarithmic) electron mass 
dependence of D through the threshold expansion technique, we could, a priori, need to 
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consider new regions to successfully obtain the subleading terms 0(m 2 /M 2 ) , 0(m 2 /(q 2 — 
AM 2 )) , etc. This is what happens, for example, if one considers the 1-loop two-point scalar 
function with one heavy mass M and one light mass m at values of q 2 > M 2 : Keeping m 
finite but smaller than any other scale present (i.e. m <C (q 2 — M 2 )/M <C y/q 2 — M 2 <C M), 
the integration region where £ 2 ~ m 2 gives a non-vanishing contribution proportional to 
m 2 / (q 2 — M 2 ). A new pattern of integration regimes should then be considered to make each 
integral homogeneous also in the m 2 scale. 

5 Conclusions 

The interest in the study of electron positron annihilation into heavy fermions has been 
ushered by the multiple features foreseen both in high-energy colliders and production at 
threshold. These include all-important aspects of the phenomenology like an accurate mea- 
surement of the heavy fermion masses (like r or t) and, the possibility, of exploring New 
Physics beyond the Standard Model. This goal requires the computation and implementa- 
tion of complete perturbative orders within the Standard Model. 

We have evaluated the QED two-photon box diagrams of Fig. ^contributing to <r(e + e~ — > 
ff) with massive final fermions (m e <C M), and we have provided a full analytical expression 
for the amplitude. Its contribution at the production threshold has also been studied and we 
have found that it is negligible because of the high velocity suppression. This non-relativistic 
analysis complements the one carried out in Ref. |7| and shows that the conclusions reached 
in that reference are not modified by the QED box amplitude. 

Finally we have analysed this low velocity behaviour using the strategy of regions to 
expand the Feynman integrals near threshold, confirming that such expansion can also be 
applied to diagrams involving heavy and light fermions. This feature allows to identify 
and evaluate, at a fixed order in the heavy fermion velocity, contributions to heavy fermion 
production or annihilation diagrams triggered by light fermions. 
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Appendix A : Integrals in the box amplitude 



In this Appendix we outline several features of the integration procedure, followed to 
evaluate the QED box diagrams, and we collect the explicit expressions for the relevant 
scalar integrals that appear in our results. 

The general structure of the two-photon box amplitude in Fig. [TJ(a), M. a takes the form 
M. a = a D + a^D^ + a^D^, where a , a M , a M „ contain Dirac algebra 7's and spinors, and 
Dq, Dp V are the integrals over the loop momentum i : 

r d 4 P 1- P ■ P P 
n . n . n _ / _ V' ZHQi (\ -\\ 

°' M ' ^ J i(2vr) 4 {£ 2 - m 2 )[{£ + k - p) 2 - M 2 ][{£ - p) 2 - \ 2 ][{P + p') 2 - X 2 ] ' 1 

which depend on three independent four-vectors and where +i5 prescriptions are understood 
in the propagators . Let us define our basis as P = p — k, Q = p + p' and R = k — k', with 
scalar products 

P 2 = t , Q 2 = s , R 2 = AM 2 - s , 

p.Q = , P ■ R = m 2 — M 2 — t , Q-R = . 

The integrals in Eq. ( |A.1| ) are invariant under the interchange {p; k} <-> {— p 1 ; —k'}. Under 
the same transformation P — ► P, Q — > —Q and R — »• R, and thus the tensor integrals D^, 
do not contain terms linear in Q, justifying our choice of basis. Tensor decomposition of 
D^ u then reads 

D li = D P P ll + D R R lt (A.2) 
£V = d pp P n p v + D PR (P^R V + R^P U ) + Drr R^R v + D QQ Q^Q u + s D 00 . (A. 3) 

Further reduction of the coefficient functions appearing in Eqs. ( |A.2| , |A.3 [) has been performed 
with the help of FeynCalc []12" |. These coefficients are thus expressed as a linear combination 



of a set of scalar integrals: D , C s , C t , Cm, B s , B t and B>m, with four (-Do), three (C a 
a — s,t, M) and two (B a , a — s,t) propagators that we collect next. 

The relevant scalar integrals have been evaluated following the method described in |13 | 



except for the rather cumbersome 4-point function D . In the latter case we have first 
calculated its imaginary part in the s-channel, following the optical theorem, and then the 
real part has been reconstructed through the t-fixed unsubtracted dispersion relation that 
satisfies D : 

^ ^ / x 1 f°° , kn.D (x,t) ,. , 

ReD {s,t) = - 4- dx , (A.4) 

7T J4A 2 X — S 

where the Principal Value of the integral is understood. We have performed its calculation 
in the K< m C M limit and, therefore, we have neglected photon masses when possible. 
As emphasized in Ref. |b|], the limit A — > is not trivial for the occurrence of terms like 



A 2 / {x — 4A 2 ), which diverge for finite Aasi-> 4A 2 but vanish for A — > at fixed x 7^ 4A 2 . 
As a consequence the photon mass should be kept finite until the final stages. 

The scalar integrals that appear in the two-photon box amplitude result in Eq. (Q) through 
the Ti functions of Eqs. (p|-|T(]|) have been evaluated in the limit where A m M and 
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for the specific cases p 2 = p' 2 = m 2 , k 2 = k' 2 = M 2 , (p + p') 2 = (k + k') 2 = s, (p — k) 2 = t. 
They read : 



Dn 



d 4 



i(2vr) 4 [£ 2 - A 2 ] [{£ + p) 2 - m 2 )] [{£ + p + p') 2 - A 2 ] [{£ + k) 2 - M 2 } 



-1 



8ir 2 s (M 2 - t) mM 



, M 2 — t -s-iS 
In In ■ 



A 2 



(A.5) 



D = £>„(*->«), 



(A.6) 



C, = 



d 4 



i(2vr) 4 [£ 2 - A 2 ] {{£ + p) 2 - m 2 } [{£ + p + p') 2 - A 2 ] 



32tt 2 s 



In 2 



-s-i<J 



7T 

y 



(A.7) 



C t = 



d 4 



i(27r) 4 [£ 2 - M 2 } [(£ - k) 2 - A 2 ] [(£ + p - k) 2 - m 2 



-1 



16vr 2 (M 2 - 1) 



aJ-L 

\M 2 



' M 2 — t s 



'M 2 -f 



Li 2 ( — ) + In 2 ( I + In [ I In 



Mm 



Mm 



m 
A 2 " 



, (A.8) 



C t = C t (t-u) , 



(A.9) 



Cm — 



i(27r) 4 [£ 2 - A 2 ] [(£ + A;) 2 - M 2 ] [(£ + k + A;') 2 - A 2 ] 



f 



16tt 2 s/3 



- 2Li 2 (l - /?) + 2Li 2 



-21n/3 ln(l + /3)+ivr In 



1-/3 



1 

2 



1 + /3 



- 2 Li 2 (-/3) 



1 + 



(A.10) 



r d D £ 
J J(2n)V 



(2ir) D [£ 2 -X 2 ] [(£ + p + p') 2 - X 2 } 



- 1 / A , -s-i<5 ' 
I A + ln 2 



16tt 2 V 



(A.lf) 
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d D £ 1 



i(27r) D [£ 2 - M 2 } [(£ + p- k) 2 - m 2 } 

-1 / -t / M 2 \ M 2 . M 2 — t \ . . 

A + ln- + ln 1-— - — In— — -2 (A.12) 



B 



lQn 2 \ p, 2 \ t ) t M 2 

d D £ 1 



M 



i(2n) D [£ 2 -X 2 ] [(£+k) 2 ~M 2 ] 
-1 



16vr 2 




(A.13) 



where ^(x) is the dilogarithm function. The two-point functions have been regularized 
within dimensional regularization in D dimensions and A = 2fi D ~ 4 /(D — 4) + 7^ — ln(47r), 
with fi the renormalization scale. From the full expressions above we see that only the 
integrals Ct, C t , D and -Do are infrared divergent for vanishing photon mass (A — > 0). 
However, as remarked in the main text, the combinations sD — 2C t (or sD — 2C t ) are 
infrared finite; accordingly all the infrared divergent contribution is provided by D and Do 
m Eq. flnD that carry a In A 2 factor. 



Appendix B : Infrared divergence of the QED box dia- 
gram 

There are several well-known facts on the structure of infrared divergences in QED that 



are relevant for our discussion [15| 



- Virtual photon radiative corrections between the external legs of a divergenceless root 
diagram generate an infrared divergent contribution that follows a specific pattern in the 
perturbative expansion. Such a structure provides a factorization of the resummation 
of the divergences at all orders. 

- All the infrared divergence in virtual photon radiative corrections commented above, 
arises from the eikonal approach in the propagator of the radiating external legs. For 
spin 1/2, for example, with k the outgoing soft photon momentum of e^k) polarization 
and p the ingoing external momentum, the modification of the fermion wave function 
reads : 

photon 1 (2p-k)-e - 

u(p) ► 2 u — ~1^ u (p) = 72 — Ti — u ^ ' ( R1 ) 

p — ft — m + id k —2k ■ p + ib 
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Figure A: Diagrams contributing to the infrared divergence of the QED box diagram. The 
wavy line corresponds to a hard photon and the dashed line to a soft photon. As explained 
in the text the infrared divergence factorizes and the spinor structure is the one of the hard 
diagram (without radiative corrections). 



that, in the eikonal approximation reduces to 

soft 2v ■ E 

«(p) — > p * n ii(p), (B.2) 

photon K — Ik ■ p + tO 

neglecting, essentially, the spin of the radiating field. 

Hence to extract the infrared divergent part of the QED box diagram in Fig. [I] we need to 
implement the eikonal approximation into the amplitude M. a in Eq. @ and the crossed M.^. 
This corresponds to evaluate the four diagrams in Fig. [A]. These are built from the tree-level 
diagram for e + e~ — > ff through one photon annihilation, by attaching a soft photon between 
an ingoing and an outgoing external leg in all possible ways. Their evaluation gives : 



M 



It = v e (p)l,u e (p)^u f (k)Yv f (k> 




(B.3) 



where infrared finite terms have not been written. In fact this result has been obtained by 
integrating over the full range of momentum of a massive photon. Rigorously we should 
define the infrared contribution by imposing an upper limit on its momentum \q^\ < A, and 
would then be replaced by A 2 in the logarithm of M{^ x - In Eq. ( |B.3[ ) we have explicitly 



nr 



stated the factorization between the hard gluon exchange, on the left, and the soft photon 
exchange inside the square brackets. 
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Alternatively we can evaluate M.{^ x from our result in Eq. (||) and we obtain 



M IR 



box 



l Q 2 f (M 2 - t)L^ + 2 L+ K - M (L+ K - L-A ( M 2 — 



8n 2 s 



M 2 — t 



Mm 



(M 2 - u)L+* 



2L^ K - M{L\ K -Ll K 



M 2 -u 



ln 



M 2 -u 
Mm 



x ln 



x 2 



— s — iS 



(B.4) 



where the spinor operators Lf K have been defined in Eq. (f|). Then, using the following 
relations : f\ 



(M 2 - t) L\ K = 2 L+ K - M - Ll K 
(M 2 -u)L^ K = - 2 L 2 K - M - L 



we finally get : 



e 2 Qf 
2s 



M z -t 



-s — 2<T 



(B.5) 



(B.6) 



whose infrared logarithm coincides with our previous result in Eq. (|B.3|) , since P K u e (p) = 
u e (p) in Lf K , being k the massless electron helicity. 

We conclude that the infrared divergence of the QED box diagram satisfies the expected 
features [TE\ and hence its cancellation should take place when real soft photon radiation 
contributions, at a fixed a perturbative order, are taken into account. 



Relations between spinor operators like these can be obtained by explicit evaluation in a particular 
reference frame or transforming the operators into traces in the spinor basis, hence working with Lorentz 
invariant expressions. 
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